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Abstract 

We explore the t'wo-loop renormalization of the specific heat for an interacting disordered elec¬ 
tron system in the case of broken time reversal symmetry. Within the nonlinear sigma model 
approach 'we derive the two-loop result for the anomalous dimension which controls scaling of 
the specific heat with temperature. As an example, we elaborate the metal-insulator transition in 
d - 2 + e dimensions for the case of broken time reversal and spin rotational symmetries and in 
the presence of Coulomb interaction. In this situation scaling of the specific heat is determined 
by the anomalous dimension of the Finkel’stein operator which is the eigen operator of the renor¬ 
malization group complementary to the eigen operator corresponding to the second moment of 
the local density of states. We find that the absolute values of the anomalous dimensions of these 
operators differ beyond one-loop approximation contrary to the noninteracting case. 

Keywords: metal-insulator transitions, nonlinear sigma model, renormalization group, 
multifractality 


1. Introduction 


The phenomenon of Anderson localization 111] has been attracting a lot of interest for more 
than 50 years since its discovery (see e.g., 0 ). The most intricate situation exists m d - 2 
dimension in which, depending on the symmetry class, a noninteracting electron system can be 
fully localized, fully delocalized or undergoes the Anderson transition. The most convenient 
tool to study the metallic phase and Anderson transition in noninteracti^ electron system is the 
low energy effective action called nonlinear sigma model (NLSM) y|-l8l]. This effective theory 
describes interaction of diffusive modes on scales larger than the mean free path. In d - 2 
this interaction results in logarithmic divergences which are summed by NLSM in a much more 
convenient way than a standard diagrammatic technique (see e.g., IT floll i. In the presence of both 


time reversal and spin rotational symmetry (class AI of the Wigner-Dyson classification 111 iMl 311 1 


the two-dimensional (2D) noninteracting electron system is believed to be always localized at 
zero temperature: there is no Anderson transition 11411 . 
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Inevitably electron-electron interaction becomes important at low temperatures. First of all, 
inelastic electron-electron scatteringwith a small (compared to temperature) energy transfer de¬ 


stroys quantum phase coherence 
quantum correction to conductivity 


This leads to a temperature (T) dependence of the 

_ Additional dependence of conductivity on T appears 

due to virtual electron-electron scattering Q. Physically, this occurs via coherent scattering 
electrons off the Friedel oscillations 0. Strong in comparison with Fermi liquid temperature 
dependence exists also in such thermodynamic quantities as the specific heat and static spin sus¬ 
ceptibility (see e.g., Ind - 2 both weak localization and electron-electron con¬ 

tributions to conductivity are logarithmic in temperature and opposite in sign. This suggests that 
the 2D metal-insulator quantum phase transition is possible in the presence of electron-electron 
interaction. 

The first extension of the one-parameter scaling theory of Ref. Cl to the case of electron- 
electron interactions was performed in Ref. lEl . Although this semi-phenomenological theory 
suffered from confusion between the local and thermodynamic density of states it put forward an 
important idea of the two-parameter scaling description of the metal-insulator transition in the 
presence of electron-electron interaction. Such multi-parameter scaling description was proven 
to be correct when NLSM has been derived for the case of an interacting electron system ll2l . 
With the help of one-loop renormalization group (RG) analysis of this NLSM an interplay of 


electron-electron interaction and disorder was analyzed Il25l - l30l] . In 2D case for the symmetry 


class AI delocalization due to electron-electron interaction overcomes weak localization in the 
weak disorder regime. This yields the metallic behavior of conductivity at low temperatures S. 
This fact supports existence of 2D metal-insulator transition in the presence of electron-electron 
interaction. 

A change in resistivity from insulating to metallic behavior with increase of electron density 
was measured in Si metal-oxide-semiconductor field effect transistor ll3lll32ll . Similar behavior 
of resistivity was experimentally observed later in a variety of 2D electron systems (for review, 
see 1 ( 3 ^ 3711. Observed temperature and electron density dependence of resistivity resembles the 
expected behavior of resistivity near a 2D metal-insulator transition and was found to be in rea¬ 
sonable agreement with the predictions of the two-parameter scaling theory ll38l - l^ . However, 
recent thermodynamics and transport measurements in Si metal-oxide-semiconductor field effect 
transistor suggest that the observed strong temperature and electron concentration dependence of 
resistivity occurs in the regime of nondegenerate Fermi system and, consequently, has nothing to 
do with the metal-insulator quantum phase transition 1141 1 - 14511 . These new experimental results 
call for development of the transport theory of nondegenerate strongly interacting 2D electron 
system, on the one hand, and more detailed understanding of the 2D metal-insulator transition 
within NLSM approach, on the other hand. 

At present, there is not much known on renormalization of the Finkel’stein NLSM beyond 
one-loop approximation (the lowest order in disorder). There are only few results within two- 
loop order approximation. The renormalization of the specific heat and static spin susceptibil¬ 
ity has been studied near the Stoner instability ilH. It was demonstrated that there is no 
metal-insulator transition in 2D electron system with Coulomb interaction and with broken time 
reversal and spin rotational symmetries 14811 . On the contrary, existence of the metal-insulator 
transition was shown in 2D interactin g el ectron system with A( —> 00 flavors (the action of NLSM 
is invariant under S\J(N) rotations) 14911 . Recently, existence of multifractality in moments of 
local density of states in the presence of interactions has been established ifsollsill . This situa¬ 
tion is in sharp contrast to the knowledge on noninteracting NLSM for which beta-function and 
anomalous dimensions of RG eigen operators are known to the fifth 052145411 and fourth ll55l457[l 
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loop orders, respectively. 

In this paper, we consider the two-loop renormalization of the specific heat for an interacting 
disordered electron system. For simplicity, we assume that the time reversal symmetry is broken, 
e.g. by a weak magnetic field. This assumption allows us to avoid contributions from the Cooper 
channel. In the absence of electron-electron interaction the system under consideration belongs to 
the symmetry class A. We derive the two-loop RG equation for the Finkel’stein parameter which 
determines the temperature behavior of the specific heat. In the case of additionally broken spin 
rotational symmetry we compare the anomalous dimension of the Finkel’stein operator in the 
NLSM action with the anomalous dimension of the second moment of the local density of states. 
We find that (i) these anomalous dimensions have opposite sign and (ii) the absolute values of 
these anomalous dimensions are different beyond one-loop approximation. 

The paper is organized as follows. In Sec. |2]we introduce NLSM approach. Next we present 
details of the two-loop computations of the Finkel’stein parameter in = 2 + e dimensions 
(Sec. O. In Sec. |4]we consider the metal-insulator transition in d - 2 + e dimensions in the 
electron system with broken time reversal and spin rotational symmetries. We end the paper with 
Conclusions (Sec. |5]). Some additional details of two-loop calculations are given in 


Appendix A 


2. Formalism 


2.1. Nonlinear sigma model action 

For the case of preserved spin rotational but broken time reversal symmetries the action of 
NLSM is given as the sum of the noninteracting part. So-, and contributions arising from the 
interactions in the particle-hole singlet and triplet channels, S int (for review, see OSSlISQIl '): 


S — S rr + S I 


( 1 ) 


Here the noninteracting part is given as 

5,. = J drTriVQf, (2) 

where g = 27TvD is the total Drude conductivity (in units e^/h and including spin). The 
Finkel’stein part of the action which involves interaction is as follows 


= r drTr[l"„crjQ]TT[lZ„crjQ]+4nTz^ fdrTr,,(Q-A) 

a,n 7=0 ^ ^ 

-27rTz^J dr Tr?j A. (3) 

Here Fq = F^ and Fi=F 2 =F 3 =F, denote the interaction amplitudes in the particle-hole 
singlet and triplet channels, respectively. The parameter Zw is frequency renormalization factor 
introduced by FinkeTstein 13 ■ We use the following matrices 

aX = sgn n (/^l = (4) 

where a,/3 = 1,..., W are replica indices. Integer numbers n and m correspond to the Matsubara 
fermionic energies e„ = jiT{2n + I) and e„, = nT{2m + I). The four Pauli matrices. 


o-\ 


0-2 

3 





= b ? 


(5) 





operate in the spin space. The matrix field Q{r) acting in the replica, Matsubara, and spin spaces 
obeys the following constraints: - I and TrQ - 0. 


2.2. ‘T'-algebra and 'F-invariance 

The NLSM action ([TJ involves the matrices in the Matsubara frequency space. Formally, 
Matsubara frequencies runs from minus to plus infinity which makes matrices of infinite size. 
To perform actual calculations with such matrices we introduce an ultraviolet cutoff for the 
Matsubara frequencies. Following Ref. we introduce additional cutoff Nm < which 
separates non-trivial and trivial (beyond which the Q matrix equals A) parts of the Q matrix. At 
the end of calculations the limit Nm, NL —> oo should be taken. 


As known SEIli, rotations of the Q matrix with a matrix exp(/;f') where^ = Tja,nXVn(^o 
play an important role. Such rotations correspond to the gauge transformations in the original 
fermionic language. In the limit Nm,NL oo and NmINL —> 0, the set of rules known as F 


(6) 


algebra ObOT allows one to establish the following relations (j - 0 , 1 , 2 , 3): 


Tr/“o-^e'^2e = Tr/“o-j2 + 

Tr rje'^Qe^'^ = Tr t/Q + ^ inx“ Tr /“o-qQ ■ 




Using Eq. one can check that, provided F^ = -z^, the NLSM action is invariant under global 
rotations of the matrix Q with the matrix exp((y) (so called F invariance}^. We remind that the 
constraint Fj = —z^ corresponds to the case of Coulomb interaction ll24ll . Since the relation 
Tj = —Za> allows additional symmetry of the NLSM action, this relation remains fulfilled under 
the RG flow. 


2.3. Thermodynamic potential 

The thermodynamic potential per unit volume is determined by the NLSM action: 

Q = -^ln JD[Q\ exp5. 


(7) 


where V stands for a sample volume. At the classical level, for Q-N, the thermodynamic poten¬ 
tial is equal to Q = —TS f[A] = InzFT'^ Tr77A. The quantum corrections to the thermodynamic 
potential determine the renormalized value of the frequency renormalization parameter 16211 : 


1 g n ^ <5^[6]) 

27rTr77A5r7’ 5f[A] ' 


( 8 ) 


We note that z!^^ is responsible for the non Lermi-liquid temperature behavior of the specific heat 

S. 


3. Two-loop renormalization of 

3.1. Perturbative expansion 

Lor the perturbative treatment (in 1 fg) of the NLSM action O we shall use the square-root 
parametrization 

g = tr-H A Vl - = n. (9) 

\W KJI 
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The blocks w and w are independent matrix variables. They have the following nonzero elements 
in the Matsubara space; w„^„^ and with wi ^ 0 and n 2 < 0. It is convenient to represent 
w and vv as the linear combinations of the Pauli matrices; ^^d = 

In what follows, we use the convention; ui, « 3 , us, ■ ■ • ^0 and n 2 , « 4 , 

Expanding the NLSM action ([1]) to the second order in W, we obtain the following propaga¬ 
tors for diffusive modes (j -0,1,2, 3); 

s 

where q stands for the momentum, n \2 - n\ — n 2 and un - Em — Sn 2 - The standard propagator 
for diffuson is given as {a>„ = 2nTn): 

^ q-+ %zj\u)n\lg- ( 11 ) 


IbTrTT 


g 




, (10) 


The diffusive modes renormalized by interaction in the singlet {!D^q\iu) = 2)*(/w)) and triplet 
{!D^q\iu) - Df\ito) - Df\ioj) = D'^iico)) particle-hole channels are as follows 

[Dg(iaJn)r^ ^q~ + 8(z^j -H r,)|w„|/g, 

( 12 ) 

[D;(/w„)]-* ^q^ + 8(z^ + T,)\u,\lg. 

For the purpose of regularization in the infrared, it is convenient to add the term to the NLSM 
action; ^ 

S + ^ J drTrAQ. (13) 

This leads to the shift of the momentum squared, q^ ^ q^ + h^, in the propagators (fTTli and (fT2li . 


3.2. One-loop perturbative results 

Before going to the two-loop results we remind briefly the one-loop perturbative results for 
Zw Expanding the NLSM action to the second order in W and using Eq. (US, we find the 
one-loop perturbative result; 

d Qf') _ 5 f[A] 

'df~ ~ fV 



Here we use the following notations; jj - Ty/z^ and = f d‘'q/{2TtY. Evaluating integral over 
momentum \t\d - 2 + e dimensions and using Eq. (I8]l, we obtain 11241] ; 


Z^^Za,\l 


th"^ 


-{js + 3r,) 


(15) 


Here we introduce resistivity t — 4Q^/g, where = Siil[2(2jtY] and Sd - 2n‘^^^IT{dl2) is the 
area of the ^f-dimensional sphere. 
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3.3. Two-loop perturbative results 

The two-loop contribution to the thermodynamic potential can be written as 


Here the term. 


c(3) _ 

int 2 


Q' 


jiT 


V 




( 16 ) 


r < 


drTTl^crjWTTr„crjAW, (17) 

J =0 a,n 

appears from the expansion of the NLSM action to the third order in W. The forth order terms 
are given as 

\j=0 ) 

16z 

X |^2/l^ H-(W56 + Wrs) - (90 + 9i)(92 + 93 ) “ (90 + QiliQl + ^ 2 )] (18) 


and 


g 


o (4) _ 

int - g 


r drTrlf.cTjKW^TrlZ, 


o-jAW\ 


(19) 


J=0 a,n 


The symbol sp denotes the trace over the spin space only. Performing contraction with the help 
of the Wick theorem and Eq. (fTOl i. we obtain 


g 


r ^ 

I ^ min{w„„ Un}[iyq\iU)n) + T)^p\i(x)m)\ ^ r jDq(iUn)D^fiiu„) 

dp’‘l m,n >0 j =0 

3 

X ^ T fDp(iaj„t)D^J \ioJm), 


( 20 ) 


/=<> 


(4K_ 47:7NrV (8 

int' 


j )7 2 “‘ 


g \6/ 

+ rj(rj + r)lDp(ia>m)lDp(i(^m) 


mm{aj,„,a>n]Dqiiaj 

m+«) 

327:TNrV 


g 


3r,(rj - r,)Dp{ito,„)D'piia>,„) 
dp^‘1 m770 


487TTNrV /8 


g \g 


I 

•JP’I m,n>0 


DpiitJm) + 'D'qiio) m+«) 


D‘p{iiD,„)Dq{iu) 

m+n ), (21) 
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and 



IbnTNrV /S j r min{w„„w„)2)p+,(!w,„+„)| — 

8 \8l Jp,qm,n>0 ^ ^ 

3 

X ^ r jDp{itJ„,)D\i’{ioj,„) - —r^rf£)^(/w„)£)p(!w„,)£)p(!w,„) 

>=o ^ 


A^nTNrV 
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( 22 ) 


Next we combine the two-loop contributions (l20l) - (1221) and rewrite them in the form 
with the help of following identity: Y^^„^Q\m\\{m,n]y{m,n) - Y^m,n>ofn[y{m,n + 
m)+y{n+m, m)]. Taking the limit w„, —> 0 and evaluating the momentum integrals and Matsubara 


frequency sum, we find (see Appendix A i 


!,2e,2 


4 = 1 + + ^^(^2 + ebi) + 0(t^) 1. 


(23) 


Here we add the one-loop contribution (fTSl l and introduce the following functions of interaction 
parameters: 


bi = -js - 3rr, b 2 = (3y, -h ys)[f{7s) + 3/(r,)] h- 6yj - 2y,, (24) 


and 


bi - 


^yt 


j=o 


2f(yj) + \n(l+yj)-2- 


il+yj)\n\l+yj) 2 + 3yj 


^yj 


yj 


M-yj) 


-6yt[yt - li 2 (-rr)] ■ 


(25) 


Here f(x) - 1 - (1 -H l/;ic)ln(l + x) and lik( 2 i:) = stands for the polylogarithm 

(Jonquiere’s function). 


3.4. Anomalous dimension of Zoj 

Since operator Tr AQ determines the local single-particle density of states, the momentum 
scale h acquires renormalization 16211 . The corresponding renormalized momentum scale h' is 
defined according to 

gh\TTAQ}^g'h'^TTA\ (26) 

where g' stands for the renormalized conductivity. Within the one-loop approximation, one can 
find lie! 


h' 


1 - ^ + ln(l -H yj)) 


2e 


i=o 


(27) 
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and mils in mils I 


g 





(28) 


Also we remind the one-loop results for the renormalization of interaction parameters 0241 l28h : 


, h^t 

Js = ys + cs,i—, 

, 

7 , = 7 t + c,,i—, 


Q,i = (1 +ys)0yt + 7s), 
Q,i = (1 +yi)(ys-yt)- 


(29) 


In order to extract the anomalous dimension of Zoj from the perturbative result (l23t . we use 
the minimal subtraction scheme 0640 . Let us introduce the dimensionless resistance I - t'h'^. 
Then using Eqs. (l27T i. (l28T l. ( l29l l and (l2Tt . we express f, jj and Za> via t, y' and 


t = (/i') y,), Jj = y' Zy/f, y^, y), Zcj = zL Z;„(f, y^, y;). (30) 

To the first order in t we find 

i Cj\t 

Z, = l+fli-, Z^^, = 1-^-, (31) 

e y. e 

where ai, Cjj, and c,,i are the functions of y', y,' now. To the second order in t the function Z-^ 
becomes 


z-' = 1 + 


bit 


i db 3 j \ 

b 2 + aibx - ^ -H e{b 2 + bx ^[/(yy) + - ln(l -i- y;)]) 

j=Q i =0 


(32) 


Here ai, ^ 1 . 2 . 3 , Cjj, and qj are the functions of y', y,. The RG equations can be found from the 
standard condition that t, yj,, and Zoj are independent of h'. In this way we obtain the following 
two-loop result for the anomalous dimension of Za>'. 


dlnza, 

dy 


= 4(6 ys, yt) = -t{ys + 3y,) - r 


(ys + 3y,) ^(< 3 ( 72 ) + 2 ii 2 (-rj)) + I2y/(yr - ii 2 (-rf)) 

j=o 


where y — Inh/h' and we introduce the function (see Refs. 51]) 


(33) 


c(y) = 2 -H 


2 H- y 

y 


ihi-y) 


+ 


1 H- y 
2y^ 


ln^(l + y). 


(34) 


For a brevity, we omitted the prime and bar signs. Equation (l3^ is the main result of our paper. 
We note that the relation 




dbi 


dy. 

>=0 'j 


(35) 


holds. This guarantees the renormalizability of the theory within the two-loop approximation, 
i.e. the absence of terms in the right hand side of Eq. (l33T l which diverge in the limit e ^ 0. 
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We note that in the limit y, » 1 (Stoner instability corresponds to y, = oo), the anomalous 
dimension (l3^ becomes ^ = -3y,t-l2y}f' in agreement with the result of Belitz and Kirkpatrick 
(see Eq. (6.60c) from Ref. 1591]). This result indicates that toward Stoner instability the loop 
expansion is controlled by the small parameter y,t 1 rather than f «c 1 as in the case of 
noninteracting electrons. The two-loop result (1331) for the anomalous dimension of Zw interpolates 
between the result of Ref. ll62[| for the case when the interaction in the triplet channel is absent, 
y, = 0 and the result of Refs. 1461 El] for the case y, —» oo. 

The following remark is in order here. As a consequence of the particle number conservation, 
the quantity z + T j has no renormalization 11241] . Therefore, the renormalization of the interaction 
parameter y^ is fully determined by the anomalous dimension 


^ = “(1 +ys)(zit,ys,yt)- 


(36) 


Thus we also derived the two-loop RG equation for the singlet channel interaction parameter y^. 


4. Scaling analysis 

As example of application of our result (l3^ . we consider the case of Coulomb interaction, 
y^ = -1 and fully broken spin rotational symmetry such that the triplet channel is absence. Since 
A time reversal symmetry is also absent, this situation can be realized in the system of disordered 
interacting fermions with magnetic impurities. In notations of Ref. ifs^ . this case is referred as 
the symmetry class “MI(LR)”. The anomalous dimension (l33l l becomes 


^j,(f) - t + i^+ -g-j + 0{t’). 


(37) 


We mention that the numerical coefficient (2 + ir 16) in front of the r term is different from the 


result (3 + tt 16) found in Ref. Il62ll . We suppose that this mismatch is due to the erroneous 


treatment of terms singular in 1/(1 + y,) in Ref. Il62r (see Appendix A i. 


The RG equation for the dimensionless resistance t is known up to the two-loop order Il48l] : 


dt 

Iny 


2t^ - AAt^ + 0{t'^), 


(38) 


where the constant 


^ 16 


139 (7r2-18)2 

12 


+ y ^(3) + (16 + ^) In^ 2 - (44 - y + 7^(3)) In 2 


+ 16^-lln42-81i4(i 


1.64. 


(39) 


Here Q ^ 0.915 stands for the Catalan constant and f(3) ^1.2 denotes the Riemann zeta. As 
usual, the zero of the y6-function,y6,(f,) = 0, determines the critical point ind - 2 +e dimensions: 


4 = |(1 


Ae)- 


O(e^). 


(40) 
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This critical point separates the metallic (f < f,) and insulating (f > t,) phases. At this critical 
point the correlation/localization length diverges 


^ - h V - t. 


t r’' 


dp, 


1 ? 


dt 

t=t^. 


= --A + 0{e). 
e 


The value of the anomalous dimension 


C = UQ = ^ J (2 - 2A + n^jb) + 0{e^) 


at the critical point determines the dynamical exponent 

A 


z-d- 


■C = 2+- + —(2A 
? 4 V 


V6)- 




(41) 


(42) 


(43) 


We remind that the dynamical exponent z determines the length scales induced by energy and 
temperature, Le ~ and Lj ~ at the critical point. Interestingly, the very same 

dynamical exponent is responsible for the deviation of the specific heat c,, from the Fermi-liquid- 
type behavior, c,, ~ Jo . 

In the considered case of the symmetry class MI(LR) the anomalous dimension of the 
Finkel’stein parameter determines the renormalization of the Finkel’stein part of the NLSM 
action. The term 5 f is an example of the operator bilinear in the matrix Q. This operator is 
T invariant, local and does not involve spatial gradients, consequently, it is also invariant with 
respect to spatial rotations of the matrix Q by matrix exp((:f'). Such spatial and time dependent 
rotations of Q correspond to the gauge transformation of the original fermions 1^ . Since the 
anomalous dimension of z^ is finite in the limit e —> 0 within the two-loop approximation, the 
operator 5 f is eigen operator with respect to the RG. We emphasize that ^~{f) > 0 in the two-loop 
approximation. 

Recently, it was shown 5J.] that the polinomial in Q operators corresponding to the 
moments of the local density of states are eigen operators of the RG. The second moment of the 
local density of states is expressed in terms of the operator 7G which is bilinear in Q similar to 
the Finkel’stein term. However, contrary to S f, the operator K 2 is not invariant due to lack 
of gauge invariance in the local density of state. For the symmetry class MI(LR) the anomalous 
dimension of the operator K 2 within two-loop approximation is given as 


(44) 


At the critical point the anomalous dimension of K 2 is negative: 

( 2 \ 2 
2-2A-^ ^+G(e'). 


(45) 


In the two-loop approximation the anomalous dimension of the q-\h moment of the local density 
of states is expressed via the anomalous dimension of the second moment, A^ = {q{q — 1 )/ 2 ]A 2 . 
They are negative and nonlinear functions of q. Thus the moments of the local density of states 
demonstrate the multifractal behavior in the presence of Coulomb interaction. We mention that 
within one-loop approximation the anomalous dimensions ^2 and are the same except the sign. 
At the two-loop order they become essentially different. 
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One can construct the polynomial in O eig en operators which corresponds to the higher 
moments of the local density of states 15(1 Isill . Similarly, one can study the higher order in 
Q eigen operators which are 'T invariant. However, at present they are not known beyond the 
operators with four Q matrices OhSll . 

It is instructive to compare the scaling results obtained above for the case of Coulomb 
interaction with the results for the system of fermions with short-range singlet interaction 
(-1 < js ^ 0). The case of short-ranged singlet interaction lies in an attraction region of 
the noninteracting fixed point. In the absence of interaction the symmetry class MI(LR) is just 
the unitary Wigner-Dyson class A. The Anderson transition in the class A is described by the 
following yS-function 152 - 541 




(46) 


The critical point and correlation length exponent are given as 






,/o) ^ 




(47) 


The anomalous dimensions of the operator K 2 within the four-loop approximation is as follows 


3f(3)^4 , 




0(e), 


(48) 


In the noninteracting case the multifractal exponent for the operator K 2 becomes 


Af =^-(r) = -(2e) 


r(0),i0) 


1/2, 


3f(3) 


^ + 0 (^ 1 ^ 


(49) 


In addition, to the operator K 2 there is the other eigen operator bilinear in Q. Its anomalous 
dimension is also known up to the four-loop order 15514571] : 


8 8 


• 0(e), 


(50) 


The corresponding critical exponent is given as 


^il = ^if(Q - (2e)02 - + 0(^1^). (51) 

We note that the anomalous dimensions and are different only by sign up to the third-loop 
order. Their absolute values become different only at the forth loop order. We mention that at the 
one-loop approximation the anomalous dimensions and for noninteracting case coincide 
with the anomalous dimensions and ^2 for the case of Coulomb interaction. Therefore, one 
can expect that the eigen operator bilinear in Q with the anomalous dimension transforms 
into the Finkel’stein operator in the case of Coulomb interaction. 


5. Conclusions 

To summarize, we studied the two-loop renormalization of the FinkeTstein parameter Zoi 
which anomalous dimension controls the scaling of the frequency and the specific heat for an 

11 















interacting disordered electron system. For simplicity, we considered the case of broken time 
reversal symmetry in order to avoid additional difficulty due to the Cooper channel. Under this 
assumption we derived the two-loop RG equation for Zoj (see Eq. (l33T l) the right hand side of 
which depends on the interaction parameters in the singlet and triplet channels. Our result (l33T l 
interpolates between the result of Ref. | 1 for jt - 0 and the result of Refs. fory, oo. 

We consider the metal-insulator transition in d - 2 + e dimensions in the case of the symmetry 
class MI(LR), i.e. with broken time reversal and spin rotational symmetries and in the presence 
of Coulomb interaction. We compared the anomalous dimension of the Finkel’stein operator S f 
in the NLSM action with the anomalous dimension of the second moment of the local density 
of states. Within one-loop approximation both anomalous dimensions coincide in absolute value 
and are equal to ones in the absence of interaction. However, in the two-loop approximation their 
absolute values deviate from each other and from anomalous dimensions without interaction. 

Finally, we mention that in order to explore the metal-insulator transition in the presence of 
spin rotational symmetry results for the two-loop renormalization of the spin susceptibility and 
conductivity are need. They will be published elsewhere. 
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Appendix A. Evaluation of two-loop integrals 

At first, we combine the two-loop contributions (l20l i - (l22l) together. Next we rewrite 
them in the form 2,„>o We used the following identity: - 

2m,n>o n + m) + y{n + m, m)]. Finally, taking the limit w„, —> 0, we find 

q(2) ^ ^ 

OJmdZoji (A.l) 

m>0 


where 


6Zoj =- 


167tTz7, 


8 


-) r + yjD„{ico,.)nf(ico,.)DgiQ) 

8 / JlFq „>o ( ^ j =0 


+ ■ 
+ . 


■ yjD{ico„)D^f{ico„)Dg{iaj„) - 2yjD^^\iu„)DqiO)Dp+g(ia)„) 

■ y]co„D\P(iaj„)Dl(Q)DpF,iioJ„) - yjDpiicj„)Dl(Q)] 
3yj[-Dp(i(L)„)Df(i(L>„)Dg(0) + Dp{iUn)J2fp\iUn)Dq{i(jJn) 

+ D^g!\ia)„)DqiQ)Dp+g{ia)n) + Dp{ia)n)Dq(Q)Dp+giia>„)^^ 
^ysZoj r 

8^ Jp,q 


mo)Dg(O). 


(A.2) 
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Now we set the temperature to zero and will study the dependence of dzoj on the momentum scale 
h only. Then we find 




l6Zoj i ’ijt + Js 


8 


yj\'^Qivh' ■^iio('yf) -^loi^Ty) 27Qj[('yj) +'yy7o2i(Tr) ■^02o(®)] 


2=0 


- 3r?[-/?io(rr) - Amijt) - 4niyt) - /o°„(0)] - \. 

Here we introduced 


(A.3) 


J^Jj) = 


8 J p,<i 


(A.4) 


and 


Jq - 


I 


®«(0) 


Appendix A.l. The integrals Jq, ond 7^jq 
U sing the result 

2Qd/rT(l + e/2)r(l - ejl) 


I 




we find 


Jq - 


AAJi^^ 


A,^Tl]Y\\-el2)Y\\+el2). 


Next, with the help of the results 


I 


8 


CO 

/‘'"I 


£)2(0) = Yldh^-^Y{1 + e/2)r(l - e/2), 


, 4Q,/r-^r(l+e/2)r(l-e/2) 

®" --(l+r,W2 + e)-■ 


we obtain 


•^02o(yj^ - 


4AM 


2e 


(l+r;)e(2 + e) 


Using the following relation between integrals 


8 


oo 

/‘'"I 


/A ln(l+-)/,■) 

^ 

T; 


I 




we find 


■^110(72) = 


AAM^" ln(l + yj) 

e2 jj 


(A.5) 


(A.6) 


(A.7) 


(A.8) 


(A.9) 


(A. 10) 


(A. 11) 
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Appendix A.2. The integral 

Next we consider the integral 7 ‘jqj. With the help of the Feynman trick (see, e.g. 
we write 

A,r(l-e)/r2^ 

Toiiyj), 




where (see Eq. (A26) of Ref. 


er 2 (l - e/ 2 ) 


roi(r;) 


= n r K 

U =1 n ) u=l 


^ -l-e/2 -l-e/2 

■^3 -^12 _ 

Xi + (1 + yj)x2 + X3 


Performing the change of variables from xi, X 2 and X 3 to s and 11 : 


1 - M 


Xi = 


3C2 = 


s + 1 ’ s + 1 ’ 

and integrating over s, we find 


X3 = 


s + 1 


0 ^ u < 1, 0^s<oo 


du 


^ 2 r"(l - e/ 2 ) r 

r^er(l-e) J 

1+yj 


2 Ei(-e/ 2 ,-e, 1 - e, 1 - m). 


Here 2 E 1 (a,j3, y, z) denotes the hypergeometric function. Rewriting the integral as 

1 


ToiCyj) 


1 

r/■ J 


l+yy 


-e 


4r2(l - e/ 2 ) 


2 Ei( 1 -e/ 2 , 1,2 +e/ 2 , m) 


(2 + e)r(l-e)^ 

and using that 2 ^ 1(1 - ejl, 1,2 + e/ 2 , m) —> - ln(l - u)lu in the limit e —> 0 , we find 

r 2 (l-e/ 2 ) 


Tmiyj) = - 


' r/Fd - e) 


21 n(l+ 7 ,) 1 , 

'^ --\n\\+yj) 


Hence, we obtain 


Jmiyj) = 


AJi 


le 


yj 


2\n(l+yj) \nHl+yj) 


2e 


Appendix A.3. The integral 7qjj 


Using the Feynman trick (see, e.g. Ref. we write 

,0 , , A,r(l-e)/ 22 ^, ^ ^ 

Jouiyj)- er2(l - e/2) 


where (see Eq.(A23) of Ref. Il48l] ) 
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l*:=lX / Va=1 
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Xk - 1 
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Ref. 1164 ). 

(A. 12) 

(A. 13) 
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(A.15) 

(A. 16) 

(A. 17) 

(A.18) 

(A. 19) 
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Changing variables from xi, X 2 , X 3 to s and u (see Eq. (IA.14b ) and evaluating integral over s, we 
obtain 


So(7j) 


-u 


du[u(\ — m)] 

(1 + yj)u + 1 — M 


2 Ei(l,-e,l-e/ 2 ;l-M(l-M)). 


(A.21) 


Rewriting the integral as 


5o(ry) 


f fifuMl - a C du[u{l-u)Yl^ ,,,,, 

7 TT—- 2 Ei(l,-e, 1 -e/ 2 ; m( 1 -«))-- —-———-, (A. 22 ) 

j (1 + ypu + 1 — M fj(i + ypu + i — u 


and using that 2 ^ 1 ( 1 , -e, 1 - e/ 2 ; m (1 - «)) ^ 1 + eln[l - m (1 - u)] in the limit e ^ 0 , we find 
the following result 


Soiyj) = — 


2 ln(l + yj) 2 


Ii 2 (-r 2 )+ 


f Jj Jj 
Here we have used the identities 

li 2 (a) - li 2 (l) + Ina ln(l - a) - - li 2 (l - a), 0 ^ a ^ 1 , 

1 2 

- Ii 2 (l/a) + li 2 (l) + ln(a - 1 ) Ina - - In a = - li 2 (l — a), a > 0 . 


Hence we obtain 


Jouiyj) 
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2\n(l+yj) 2 
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Appendix A.4. The integral /qji 

As above we use the Feynman trick (see, e.g. Ref. to write 

,1 A,r(l-e)h2^,;^_ ^ 

■^02l(‘^) - „2/i 


where 


s\{yj) 


n f 

U=1 \k=l 


Xk - 1 


eFHl-e/l) 


(xiX2 + X2X3 + X3X1) * 


X2 


((1 +r;)xi +X3)2 


(A.23) 


(A.24) 


(A.25) 


(A.26) 
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Changing variables from xj, X 2 , X 3 to s and u (see Eq. (lA.Hl i') and evaluating integral over s, we 
obtain 


sliyj) = 


1 
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du [m (1 — m)]* 
e{2 + e) J ((1 + yj)u + 1 - m)" 


' 2 F\( 2 , -e, 1 - e/2, 1 - m (1 - m)). 


(A.28) 


Rewriting the integral as 
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and using that 2 F\{\, —e, -e/2,1 - m(1 - «)) —> [2 - m(1 - m)]/[m(1 - u)] in the limit e —> 0, we 
obtain the following result 


S\(7j) = 


2 2yj-(2 + yj)\n(l +yj) 


2 ^ 2{2 + yj)ln{l+yj) 


l+yj 


2(2 + yj) 


M-7j) + ^lnHl+yj) 


(A.30) 


Hence, we find 




2(2 + ry)ln(l+ri) 


yy 


(1 +rj)e 


7]e 


2(2+ yj) 


y]e 


M- 7 j)+^ln\l+yj) 


(A.31) 


Finally, substituting the results for the integrals into Eq. (IA.3I) we obtain the result 
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